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Abstract A new concept analogous to global hyperbolicity is introduced, based on test fields. It is shown that 
the space-times termed here "curve intcgrable" are globally hyperbolic in this new sense, and a plausibility 
argument is given suggesting that the result applies to shell crossing singularities. If the assumptions behind 
this last argument are valid, this provides an alternative route to the assertion that such singularities do not 
- - - violate cosmic censorship. 
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Penrose's strong cosmic censorship hypothesis [7] postulates that, subject to genericness conditions, space- 
0h! time will be globally hyperbolic: i.e. strong causality holds and J~^{p) n J~{q) is compact for all p,q. The 

context in which this is usually discussed [6] is that of space-times where the metric is C^~ (ensuring the 
existence of unique geodesies in the classical sense). As result, a Lorentz manifold with a metric that fails to 
^o ' be C^~ at one point p has to be viewed as a space-time from which p is to be deleted, which usually results 

in a failure of global hyperbolicity and thus a breakdown in cosmic censorship. 

(N ■ 

^ , On the other hand, there are increasingly many examples emerging of space-times that violate cosmic cen- 

(y^ ' sorship in this way, but where there is a well posed initial value problem for test fields. Global hyperbolicity 

CO . is sufficient, but not necessary for this. The physically meaningful condition is not global hyperbolicity, but 

^^ ' the well-posedness of the field equations. This suggests that we need to redefine the notion of hyperbolicity 
(or equivalently, of cosmic censorship) so as to make direct reference to test fields, in situations where a 

r^^ ' putative singularity p is an internal point of a Lorentz manifold (M, g) with g not being C^~ at p. (A 

Q^ , general discussion of this idea is given in [2].) 

O . We shall suppose throughout that we are working in a region in which there exist local coordinates (we 

^ \ regard M C IR ) in which gij and g^^ are bounded but not necessarily continuous. This is not essential for 

5-H . many of the points discussed below, but it is not my aim here to discuss global issues. The definitions that 

Zv', follow below are then made with respect to a choice of some particularl type of physical field 4> (e.g. massless 

scalar) satisfying L(f) — for a second order differential operator L — L{g). 
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There is an underlying conflict (which will not be resolved here) between the view of geometrical general 
relativity, in which one deals with differentiablity in the four-dimensional manifold, and the 3-1-1 setting 
appropriate to the analysis of hyperbolic equations, in which function-spaces are defined on three-dimensional 
hypersurfaces. The ideal would be an integrated approach in which the full coupled Einstein-matter equations 
were handled in a manner consistent with their hyperbolic nature. In this paper I am the inhabiting a helf- 
way house, in which the metric is being is being handled in terms of 4-D differentiability, while the matter 
is being regarded as a test field used to probe the metric and is described in 3 -I- 1 terms. Thus for the 
case being considered where L is second order, we shall take a foliation of the region which establishes a 
particular diffeomorphism with IR x IR and regard the field (f> as a map $ : 1 1-^ (</)(., t), (/)(., t)) taking values 
in an appropriate function space (defined by the norms in the next section) on IR , and L will have the 
form d^/dt + A^ for a three-dimensional differential operator A. For simplicity of notation, however, I shall 
usually not distinguish explicitly between $ and 0. 

The extension of the "usual" definition of an operator to the case of a non-smooth metric is not always 
unambiguous, or even possible. In the case, however, of the wave operator which will be treated here we can 
regard <I>(i) as lying in H^(Si ) x iJ°(IR3) (iJ* being the Hilbert space of i times differentiable functions. 
For smooth g and for V', e C^{JR'^) x C^(IR'*) we have that 



The IR integral on the right hand side defines, for fixed (/), a linear operator on ip which is bounded on 
-ff^(IR ) X if°(IR3), and is well defined for a general bounded invertible g. Hence there is an element A{^{t)) 
of this space such that 

thus defining A (an unbounded operator on a dense domain) for general g. 

In this context, I shall call cj) a "solution" to Lip = if there exists a foliation with respect to which 
d^/dt + A^ = 0. Note that this therefore does not imply that </) is C^. A solution in this sense is also a 
weak solution in the sense that it has a locally integrable weak derivative (a distributional derivative that is 
a function) (f>^k satisfying 

JV^9"'<t>,kX.dV^0 (1) 

for all test functions x- 

I then make the following defintions. 

M is L-globally hyperbolic if there is a spacelike surface S such that there is a 1-1 correspondence (defined 
by taking the foliation for $ to include S) between Cauchy data on S (satisfying only local constraints, if 
any) and global solutions to Lcf) = 0. 

A point p in M is L-nakedly singular ii it has no L-globally hyperbolic neighbourhood. Otherwise it is called 
L-inessential. 

For simplicity of the later exposition I shall take Cauchy data that is C^ x C^. This restriction on the 
initial conditions is stronger than is really required: the aim is to illustrate principles, not to obtain the best 
possible result. 

2. Curve integrable space-times 

Perhaps the most interesting of the singularities where the differentiability falls below C^" are the shell- 
crossing singularities, which are not too unrealistic physically and may be tractable analytically. With a 
view to showing that these are inessential, I shall prove the following: 

Theorem. Suppose given (M, 17) and a point p in M such that 

(a) gij and g*-' arc continuous 

(b) ffy isC'^ inM\J+{p) 

(c) weak derivatives gij,k exist and are square integrable on M 

(d) the distributionaUy dcRned W jki is a function 

(e) there is a non-empty open set C C IR , and positive functions 

M, N : IR^ -^ IR^ such that, if 7 is a curve with d-y/ds e C for ah s then 

(i) 7 is future timelike 
(ii) the integrals 

/•a /•a 

I^{a):^j \T]k{l{s))\^ ds and J^{a) -.^ iR'^ailis^ds 

JQ Jo 

are convergent, with 

Iy{a) < M{a), J^{a) < N{a) (2) 

and M{a) -^ 0, N{a) ^ as a ^ 0. 
Then p is \I\- inessential (where D is the wave operator). 



The conditions (a) and (c), introduced by Geroch and Traschen [5], are the minimal conditions for R^ jki to 
be definable as a distribution by the usual coordinate formula in terms of gij and g"^^ . The set C defines 
a range of timelike directions which are transverse to any shocks or caustics that may be present (this is 
illustrated in the next section). We refer to (e) by saying that the space-time is curve- integrable. 

I conjecture that the same result holds if the definition of / is altered to involve the simple modulus of F, rather 
than the modulus squared. The proof here, however, requires the stronger condition above, which essentially 
asserts that the quadratic and the linear (in F) parts of the Riemann tensor are separately integrable along 
the curves considered here. This form of the condition implies, of course, the square integrability stated 
separately in (c). 



Proof 

We will be able to apply standard theorems on the existence and uniqueness of solutions to linear partial 
differential equations, provided that we can establish an energy inequality for the solution. This in turn will 
require us to find a vector field t whose covariant derivative is known to be bounded on spacelike surfaces. 
We achieve this by taking for the field the tangent vector to a suitable congruence of geodesies. The necessary 
steps are, therefore, to show the existence of the congruence and then to compute the covariant derivative 
of its tangent vector. 

2. 1 There exists a congruence of timelike geodesies whose tangent vector has an essentiaUy bounded covariant 
derivative 

Proof of 2.1 

Since the definition of 'inessential' is local, we can shrink M to a smaller neighbourhood of p if necessary. 
Since the metric is continuous we can take a rotation of coordinates so that the surfaces St = {x \ x'^ — t} 
are spacelike in M (shrunk if necessary). Suppose moreover that p is at i = 0. Choose a fixed vector tg e C 
and let S be the spacelikc surface St^ , where ti < is to be determined later. We denote the coordinates x" 
(a — 1, 2, 3) on S* by 2/, and let tj;o denote the vector with components tg at the point on S with coordinates 

Next, we need to establish that the conditions on the metric are sufficient to ensure the existence of geodesies 
with initial tangents tijo- we shall take these to form the required congruence. Let S* : IR — > IR be 
a C°° smoothing function, i.e. S{x) > 0, S{x) — S{—x), J S{x)dx = 1, support(S') compact. Define 



and let fJ^^' = S„ 



* ^jk (where * denotes convolution). For each j/ e IR 



let K^ be the 



Sn{x) — n*S{nx) 

p(")_geodesic with initial tangent vector igg and let tl"'(s) be its tangent vector at parameter s. 

We note that, for small enough s, say s < si, these tangent vectors have components in C (with si being a 
uniform bound, independent of n). This follows from the geodesic equation 



together with 



ds 



= -F 



{n)i,{n)j.{n)k 



jk 
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(3) 



r%^'iKs'))ds' 



Sn{z) 



r)^{x{s') + Z)ds' 



dz 



< Misf's^''^ 



(4) 



for any curve with A e C, which establishes a uniform bound for the right hand side of (3). Indeed, if 

r = mindist(tl"^(0),C") 



.(«) 



(where C" is the complement of C) and q = max^; |tj-" (0)| then it suffices to take si such that M(si)^"s^" < 

r/{8{q + r)^) which will ensure that |tl" — tj;o| < r. The time ti can now be specified as sufficiently close 
to to ensure that p is covered by the curves up to si. 



We now examine the connecting vector Y, the basic tool being the geodesic deviation equation for the 
variation with y oi k" . Let (e- )a=o...3 be a parallely propagated co-frame on k- coinciding with the 

(n) 

coordinate basis at s = 0, with e-j, the corresponding frame, and define 






The geodesic deviation equation is then 
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ds 
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ds 
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Condition (e) imphes that e^ i are uniformly bounded, as are t~ , so that (5) gives 



^2yHa 
ya 



ds^ 



— ^^ II ya I 



for some constant Q, where 



a:= sup |i?("^-fez| 

i,j,k,l 



It follows that \\Y\\ is bounded by the solution z of the majorizing equation 

d^z 



ds^ 



Kaz 



subject to z(0) = 1, dz{0)/ds ^ ||r ;;^-(ti,y)|| =: V (defining V) 



in)., 



(5) 



(6) 



This equation will imply that z, and hence ||y||, can be bounded in terms of the integral of R^"'. Now the 
significance of the conditions (e) is that this integral can be bounded in terms of the integral of R. Indeed, 
since the integrals of the linear and quadratic parts of R are separately bounded, we have inequalities of the 
following form (with a, 6, ai, a2 constants) 

|a.r5^°Ms'< /Vn*9.r;.^^°Ms'< f dzSn{z) r\d,r%i^if{s') + z)\ds'<a{j{s) + i{s)) 



\Tf^'\\Ttl'\ds' < j dzSn{z) j dySniy) j^ rf^'|r;-fe(4"^(5') + ^)l|rLp(4"'(*') + y)l 

< hl{s) 
|i?^"^feiMs' < / ads' < aiM(s) + a2N(s) =: Mi(s), 



-^,<2Qa 



and hence 



say. To estimate the solution to (6) we then note that if S2 is the first value of s at which |z| = 2 (possibly 
S2 = oo) then before s^ 

fz 

ds^ 



leading, for < s < S2, to 



Thus if we choose sq so smah that 



z<l + Vs + 2Q I Mi{s)ds. 

/So 
Mi{s)ds < 1 



then we wih have z < 2 up to sq, and hence H^^^ || < 2 in this intervaL 

This bound iniphes that the function k : {y, s) i— > k"{s) is equicontinuous, and so there is by Arzela-Ascoh a 
subsequence that tends to a hniit. Choosing this subsequence gives meaning to the idea of hmiting geodesies. 

Having estabhshed this, essential boundedness of the derivative of the tangent vector follows in a similar 
way. If X' = X"ri"^' + XH''"'^' then 

y !j y ya jJ ja\'y/ U / y ■> y y yet 



We can now take the limit of this in L°^ to obtain essential boundedness. 

2.2 Solutions of the wave equation satisfy an energy inequality (12) 

Proof of 2.2 

The technique closely follows the account of Hawking and Ellis [6] , section 7.4 

We are concerned with solutions (cf (1)) to the wave equation 



(7) 



Suppose initially that (j) is C^ and define 

We let C/ be a compact set bounded to the past by S (i.e. I^ (U) n dU C S) and to the future by a spacelike 
surface i/ = /+(t/)n at/. 

Working locally, from the continuity of g^^ we can choose U to be foliated by C°° spacelike surfaces S^ . (We 
take < T < Ti.) 

Set Ur = Ur'<r ^r' and define 



E{t):= / S'HiUjy/^d' 



sy 



where n is the future normal to S^ . Our aim is to estimate the norm 



1 

\S,T 






1/2 



which is related to E by 



for positive k, K. We also introduce 



kE{T)<{U\\l,f<KE{i 



(8) 



< 



Cy]{<^,^? + 4>^)d^x 


mit'fdt' 


1/2 



1/2 



(9) 



and 



Stokes' theorem yields 



By direct calculation 



(f>'')d^X 



U-r 



1/2 



U^ 



{S'H,),jy/^d'^x -(-/+/ )S'H,njy/^d?x. 



S JH. 



(10) 



and so the left hand side of (10) becomes (with g^^(t>:ki = 0(f)) 

Estimating all the terms by the bounds available gives 

1 dl^llU' < E{r) < EiO) + cmlWL^'l + c'imlf 

for constants c, c'. If D^ — weakly, this becomes 

1 



K 



while if E{0) = 0, (11) and (9) give 



{mi,y<Eit)<Eio)+c'm\i) 



(11) 

(12) 
(13) 



E{r)<c,{\\m'^r 

for some constant ci. 

2.3 There exist unique solutions to the wave equation for C^ x C^ initial conditions 

Proof of 2.3 

We briefly recall the standard arguments (see, for example [4]) which allow us to deduce uniqueness and 
existence of solutions from an energy estimate, using the symmetry of the wave operator. 

Uniqueness is immediate: if the difference between two solutions is zero on 5* then E{0) = and (13) then 
implies that the solutions are (pp) identical. 

Let Vi be the subset of L'^{U) consisting of C°° functions that are zero with a zero normal derivative on 
H := Sri s-iid let Vq be the subset of L'^{U) consisting of C°° functions that are zero with a zero normal 
derivative on S := Sq. Then the same uniqueness argument holds for both these data conditions, and we 
have equation (13), which implies there exists a constant C2 such that 



<c2\\m 



(14) 



for G T^ (from now on all norms are in L'^). 

To prove existence subject to conditions (p = (f)Q,(j) = (f)i on S, choose a C^ function / satisfying these 
conditions and look for a function ijj ^ (j) — f satisfying zero boundary conditions and Lip = —Lf ~: x- The 
required fucntion ip will satisfy 



ipLwd X — I x'^d X 
u Ju 



for all w ieVi. From (14) 



Xwd'^x 



<C2||xlll|i^ll 
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Figure 1: space-time diagram of world lines of shell-crossing dust 



so that the map k : Lw 



Ju X'^d,'^^ is a bounded linear functional on LVi. But LVi is dense in L^, because 

■741, 



Vb is dense and if /„ (pLwd x — for all Lw e iVi and (j) £ Vq then we must have cf) — 0. So k defines an 
element xp of L^ such that J„ xpLwd^x =■ k{Lw) = Jjj xwd*x; so that cf) ^ f + ip is the required solution. 

This concludes the proof. 



3. Application to dust caustic (shell crossing) space-times 

Though there is as yet no rigorous proof, there are very strong indications [3] that shell-crossing spherically 
symmetric dust configurations produce relativistic solutions in which the flow lines of matter produce a 
caustic, as in the gravity-free case (see figure 1). 

In what follows I shall be assuming the existence of such a space-time, in which the general form of the 
matter density is the same as that in the gravity free case. Since the matter density determines the Riemann 
tensor and the connection via simple integrals in this case, we can pass from the density to the Riemann 
tensor immediately. 

From catastrophe theory, the generic gravity free caustic is diffeomorphic to the following canonical form. 
If (r, t) are the essential coordinates in a spherically symmetric situation, and the flow lines of matter are 
parametrised by t, with v :~ dr/dt constant on each line, then the lines when lifted to curves 1 1— > (r, t, v) in 
IR (thought of as a reduction of the tangent bundle) rule the surface E with equation r = f{v,t) :— vt — av^ 
for a constant a. The projection in the tangent bundle corresponds to the projection p : (r,t,v) <—> {r,t), 
and the critical points of p|S constitute the curve t = iav^ in E. The caustic is the image under p of 
this critical point set, namely t = 3a^'^{r/2)'^'^. Each point on the ruled surface makes a contribution 
P = Po{v){df /dv)^^ = po/{t — iav^) to the total density at the corresponding point of space-time, for some 
function po{v) giving the density distribution in velocity-space. 

Integrability of the Riemann tensor along curves will depend on its behaviour near the caustic. Consider, 
therefore, a coordinate straight line cutting the caustic at the image of a point on the critical point set in S 
with velocity uq, i.e. at tq — 2afQ, ig = Sawg, the line being x — xq + X{t — to). Setting t = Iq + t, v ^ vq + i/ 
and working to lowest significant order in r and i^, we obtain 

\ \ 1/2 

3avo J 
for Wo 7^ (where the condition fo > A is required for transversality to the caustic) and 



At 
a 



1/3 



7 



for Wo = 0. The key point arising from this as a consequence is that p is integrable along the curve, a result 
which is diffeomorphism invariant and so applies to the generic caustic. 

Passing to the relativistic case, as previously noted we assume that this behaviour of the density still holds, 
specifically when the metric is presented in coordinates linearly related to double null coordinates. (As 
described in detail in [3] , the choice of coordinates becomes significant in general relativity, as opposed to 
the Newtonian case, because coordinate transformations - for instance, from curvature coordinates to double 
null coordinates - are typically specified by geometrical conditions and so are not, in this case, C°°.) With 
this assumption, the Riemann tensor is curve- integrable and (see the treatment of these coordinates in [1]) 
the connection coefficients are bounded. It would then be the case that the caustic is not a D-essential 
singularity, so that cosmic censorship is not violated. 
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